function plotCSregressions(momentsModel,momentsData)
FontSizeValue = 12;
lineWidthValue = 2.5;
markerSizeValue = 10;

h_index = [1,3,6,12];
for i=1:length(h_index)
    name = ['h',num2str(h_index(1,i))];
    figure('Name',['Campbell-Shiller regressions: h = ',num2str(h_index(1,i))],'NumberTitle','off')
    % Intercepts coefficients
    subplot(2,2,1)
    hold on
    plot(momentsData.(name).maturities/12,momentsData.(name).LPYI.CSBetta(1,:),'-k')
    plot(momentsModel.(name).maturities/12,momentsModel.(name).LPYI.CSBetta(1,:),'-xk','MarkerSize',markerSizeValue)
    plot(momentsModel.(name).maturities/12,zeros(length(momentsModel.(name).maturities),1),'--k')
    xlabel('Maturity in years','FontSize',FontSizeValue)
    set(gca,'FontSize',FontSizeValue)
    hold off
    title('Intercepts')
    hold off
    
    % Intercepts: regime-dependent
    subplot(2,2,2)
    hold on
    plot(momentsData.(name).maturities/12,momentsData.(name).LPYI_Regime.CSBetta_exp(1,:),'-b')
    plot(momentsModel.(name).maturities/12,momentsModel.(name).LPYI_Regime.CSBetta_exp(1,:),'-xb','MarkerSize',markerSizeValue)
    plot(momentsData.(name).maturities/12,momentsData.(name).LPYI_Regime.CSBetta_rec(1,:),'-r','LineWidth',lineWidthValue)
    plot(momentsModel.(name).maturities/12,momentsModel.(name).LPYI_Regime.CSBetta_rec(1,:),'-xr','LineWidth',lineWidthValue,'MarkerSize',markerSizeValue)
    plot(momentsModel.(name).maturities/12,zeros(length(momentsModel.(name).maturities),1),'--k')
    xlabel('Maturity in years','FontSize',FontSizeValue)
    set(gca,'FontSize',FontSizeValue)
    hold off
    title('Intercepts: Regime-dependent')
    hold off
    
    
    % Slope coefficients
    subplot(2,2,3)
    hold on
    plot(momentsData.(name).maturities/12,momentsData.(name).LPYI.CSBetta(2,:),'-k')
    plot(momentsModel.(name).maturities/12,momentsModel.(name).LPYI.CSBetta(2,:),'-xk','MarkerSize',markerSizeValue)
    plot(momentsModel.(name).maturities/12,zeros(length(momentsModel.(name).maturities),1),'--k')
    xlabel('Maturity in years','FontSize',FontSizeValue)
    set(gca,'FontSize',FontSizeValue)
    legend({'Data','Model'},'Orientation','horizontal',...
        'Position',[0.234558531746033 0.959920634920635 0.117291666666667 0.0331632653061225],...
        'FontSize',FontSizeValue)
    hold off
    title('Slopes')
    hold off
    
    % Slope coefficients: regime-dependent
    subplot(2,2,4)
    hold on
    plot(momentsData.(name).maturities/12,momentsData.(name).LPYI_Regime.CSBetta_exp(2,:),'-b')
    plot(momentsModel.(name).maturities/12,momentsModel.(name).LPYI_Regime.CSBetta_exp(2,:),'-xb','MarkerSize',markerSizeValue)
    plot(momentsData.(name).maturities/12,momentsData.(name).LPYI_Regime.CSBetta_rec(2,:),'-r','LineWidth',lineWidthValue)
    plot(momentsModel.(name).maturities/12,momentsModel.(name).LPYI_Regime.CSBetta_rec(2,:),'-xr','LineWidth',lineWidthValue,'MarkerSize',markerSizeValue)
    plot(momentsModel.(name).maturities/12,zeros(length(momentsModel.(name).maturities),1),'--k')
    xlabel('Maturity in years','FontSize',FontSizeValue)
    set(gca,'FontSize',FontSizeValue)
    legend({'Data: Boom','Model: Boom','Data: Recession','Model: Recession'},'Orientation','horizontal',...
        'Position',[0.544454365079364 0.96374716553288 0.38625 0.0331632653061225],...
        'FontSize',FontSizeValue)
    hold off
    title('Slopes: Regime-dependent')
    hold off
    
    
    
end



end